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Abstract 

The influence of the small scale "cellular" structure of premixed flames on their evolution at 
larger scales is investigated. A procedure of the space-time averaging of the flow variables over 
flame cells is introduced. It is proved that to the leading order in the flame front thickness, the form 
of dynamical equations for the averaged gas velocity and pressure, as well as of jump conditions for 
these quantities at the flame front, is the same as in the case of a zero-thickness flame propagating in 
an ideal fluid at constant velocity with respect to the fuel, equal to the adiabatic velocity of a plane 
flame times a factor describing increase of the flame front length due to the local front wrinkling. 
As an application, the large scale evolution of a flame in the gravitational field is investigated. 
A weakly nonlinear non-stationary equation for the averaged flame front position is derived. It 
is found that the leading nonlinear gravitational effects stabilize the flame propagating in the 
direction of the field. The resulting stationary flame configurations are determined analytically. 
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I. INTRODUCTION 



Propagation of plane flames in gaseous mixtures is well known to be unstable. An efficient 
way of investigating this instability is to consider the flame front as a surface of disconti- 
nuity, expanding all quantities of interest in powers of Lf/X, where Lf is the flame front 
thickness, and A characteristic length scale of a flame perturbation. The leading term of the 
perturbation growth rate expansion has the form 

o- = coy (1) 

where Uf denotes an adiabatic velocity of a plane flame front with respect to the fuel, and 
c o — c o($) a function of the gas expansion coefficient 9 defined as the ratio of the fuel density 
(p u ) and the density of burnt matter (pb), 6 = Pu/Pb > 1- According to Refs. 

, CO(0) 

has positive values for all 6, implying an unconditional instability of zero-thickness flames, 
the Landau-Darrieus (LD) instability. In the next order in Lf/X, account of the transport 
processes inside the flame front modifies Eq. ([!} to 

where ci depends on 9 as well as on the ratio of the heat and mass diffusivities (the Lewis 
nu m be r ) S fl 3. The product c t L, , K the S o-ea„ed cut-o ff wavele ngt h, is the 8 ho r t 
wavelength limit of unstable perturbations. By the order of magnitude, A c represents also 
the characteristic length of the so-called cellular structure of the flame front, which is formed 
eventually as a result of the nonlinear flame stabilization. For many flames of practical 
interest, c\ = 15 — 20. It is the fact that Lf is relatively small in comparison with A c which 
underlies the above point of view on flame dynamics. 

Consider an arbitrary initially smooth front configuration. As a result of the rapid growth 
of the unstable flame perturbations with wavelengths ~ A c , the flame front becomes corru- 
gated within the time interval 

A t ~|. (3) 

Dynamics of the short wavelength modes are mainly determined by the transport processes 
inside the flame front, and are affected only slightly by the large scale flow. One can say 
that the small scale cellular structure of the flame front develops on the "background" of its 
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front position 



FIG. 1: Schematic representation of the scale separation. 

smooth large scale configuration (see Fig. 1). It follows all developments of the background, 
Eq. (jSJ) playing the role of the characteristic time of cell adaptation to the large scale front 
evolution. 

In practice, it is the large scale evolution of the flame front, rather than its exact local 
structure, which is often of the main concern. In this respect, an important question arises 
about the reverse influence of the flame cellular structure on the front evolution at scales 
much larger than A c . More precisely, one can state the problem as follows. Imagine that we 
have smeared the small scale rapid variations of all relevant quantities by averaging them 
over many flame cells. Then the question is what equations governing dynamics of the 
averages are. 

In connection with the above statement of the problem it should be noted that the exact 
cellular structure of flames is actually unknown. This is because the process of cell formation 
is essentially nonlinear, in the sense that it cannot be treated perturbatively in principle, 
which is the main reason of lack of its theoretical description. Only in the case 9 — > 1, 
which i8 practical , can this s_ be determined ana^ flflflflQ. 

The question of principle, therefore, is to what extent the large scale dynamics of averages 
depend on the exact local flame structure in the regime of fully developed LD-instability. 

The main purpose of the present paper is to show that to the leading order in the ratio of 
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A c to the characteristic length of the problem (Lq), dynamics of the averaged quantities are 
actually independent of particularities of the local flame structure. The latter determines 
essentially only one parameter characterizing the large scale evolution - the effective normal 
velocity of the flame front. 

Perhaps, it is worth to explain the essence of the problem in a little bit more detail. The 
above point of view on the flame propagation is based on the possibility to separate the 
local cellular dynamics from the large scale evolution of the background. This possibility is 
underlined by the following common property of the transport processes. From the mathe- 
matical point of view, all these processes are of higher differential order than those governing 
dynamics of an ideal fluid. Therefore, their relative role increases at smaller scales. In partic- 
ular, in the limit Lf — > 0, cell formation is completely determined by the transport processes 
inside the flame front. On the contrary, the role of these processes at scales L ^> Lf is 
relatively small. It should be fully realized, however, that this reasoning is inherently linear. 
It tacitly assumes that if every quantity of interest, say A, is represented as a sum of its 
averaged value (A) = Aq and the small scale fluctuation A%, then the dynamics of A^s can 
be determined solely in terms of Ao's themselves. Because of the high nonlinearity of basic 
equations governing the flame propagation, this assumption is far from being self-evident. 
For instance, averaging of a cubic combination of A's gives rise to a term (Af)Ao comparable 
with Aq, since the small and large scale parts of the flow variables are generally of the same 
order of magnitude. Clearly, equations for Ao's involving such terms would not be of great 
value, since the local flame dynamics, and therefore, the coefficients (Af), are unknown. The 
main result of the present work is the proof that such terms actually do not arise in the 
leading order with respect to A c /L . One can say that the governing equations for the quan- 
tities Aq and A\ decouple from each other. The proof consists of two parts corresponding to 
decoupling of the flow equations in the bulk, and decoupling of the jump conditions at the 
flame front, given in Sees. Ill Bl and IIIC| respectively. As an application of this result, the 
problem of nonlinear front stabilization in a gravitational field will be considered in Sec. IIHI 
The results obtained are discussed in Sec. IIVI 
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II. THE DECOUPLING THEOREM 



A. The averaging procedure 

Let us begin with the precise formulation of the averaging procedure. Denote L the 
characteristic length of the problem in question. For instance, L can be the tube width, in 
the case of a flame propagating in a tube, or be related to an external field acting on the 
system. In practice, this length largely exceeds the flame cell size, 

Assuming this, let us choose a length L satisfying 

A c < L < L . 

Analogously, denoting the characteristic time interval by T , and noting that 

rri Lq 

J n ~ — , 

Uf 



we can choose T ~ L/Uf such that 



-f < T < T . 



Given a function A(x,t), we define its space-time average over {x,t : x e (x ,x + Ax), 
Ax t = L, te (t ,t + T)} 

( A ) = JJfJ dt J d " x A ( x ^) = M*o,t ) ■ (4) 

T V 

By the definition, (A) varies noticeably over space distances ||Acco|| ~ L , and time intervals 
At ~ T . The function A thus turns out to be decomposed into two parts corresponding to 
the two scales, L and A c : 

A = A + A 1 , (A 1 ) = 0. 

As was mentioned in the Introduction, flame dynamics can be analyzed in the framework of 
the power expansion with respect to the small ratio Lf/L = e (or equivalently, with respect 
to A c /L , since A c = 0(L{)). Thus, we write A and A± as follows 

A = A^+eA^ + ..., A l = Af ) +eA^ + ..., 



dots denoting terms of higher order in e. In this notation, the large scale flame dynamics in 
zero order approximation with respect to e are described by the quantities A^\ Our main 
purpose below will be to investigate coupling between and Af\ A± , etc., and to obtain 
effective equations governing dynamics of A^\ Accordingly, all quantities will be measured 
in units relevant to the large scale dynamics. Namely, space coordinates x and time t are 
assumed to be normalized on L and L /U{, respectively. Furthermore, Uf will be taken as 
the unit of gas velocity v, while p u Uf as the unit of gas pressure p . For future reference, let 
us write down their expansions explicitly 

V = V + «! , p = p + pi , (5) 

(0) , (1) , (0) , (1) , ta\ 

v = v y ' + ev y ' -\ , v 1 = v{ + ev\ H , (6) 

(0) i (!) i (0) , (1) , fn\ 

Po = Po + m ' + ■ ■ ■ , Pi = p\ + ep\ + • • • . (7) 
Within our choice of units, we have the following estimates 

v = Oil), «i = 0(l), Po = 0(l), Pl = 0(l), (8) 

/f- = 0(l), /^- = ofi), i,t,/=l,2,3. (11) 

Let us now proceed to the examination of flame dynamics in terms of Wqi jPoi ■ 

The main result concerning the large scale flame dynamics which will be proved below 
can be expressed in the form of the following 

Decoupling theorem: The large scale dynamics of a flame are unaffected by its local cellular 
structure up to a rescaling. More precisely, the form of dynamical equations for , p^ , as 
well as of jump conditions for these quantities at the flame front, is the same as in the case 
of zero-thickness flame propagating in an ideal fluid at constant speed it Uf with respect to 
the fuel, the number il > 1 describing the flame front length increase due to the local front 
wrinkling. 

The proof consists of two parts corresponding to decoupling of the flow equations in the 
bulk, and decoupling of the jump conditions at the flame front, presented in Sees. Ill Bl and 
III CI respectively. 
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B. Decoupling of dynamical equations 



For definiteness, we will assume in what follows that external field acting on the system 
is the gravitational field, denoting its strength by g. Accordingly, L Q will be identified with 
the characteristic length associated with this field: 1 



L 



El 
\\g\\ 



(12) 



Then the dimensionless velocity and pressure fields obey the following equations in the bulk 



dv 
~dt 



where 



div v 
(vV)v 

G = 



Q 



-Vp + G + ePrAv. 



(13) 
(14) 



\G\ 



q is the fluid density normalized on the fuel density p u , and Pr the Prandtl number repre- 
senting the ratio of viscous and thermal diffusivities, Pr = vj\- 

Substituting expansions (0-0 into Eqs. (|T3*|) . (jlljl . taking into account the estimates 
(|Sjl- (fH)l . and extracting 0(1/ e) terms yields 



div v\ 



(o) 







(0) 



Of 



v, 



(0) 



+ v 



(0) 



(0) 



-VpS 0) + ePrAvr ■ 
Q 



,(°) 



Next, collecting 0(1) terms gives 



(o) 



+ e—^— + 



Of 



Of 



v, 



(0) 



+ v 



(0) 



div v, 



V\v® + e 



(0) 



+ ^! 0) ] V) v[ 



(i) 



g 



(0) 



-Vp\ 

g 



(i) 



G + e 2 PrAv\ 



(i) 



(15) 
(16) 

(17) 

(18) 



Equation (j!8)l involves both slowly and rapidly varying terms. The slowly varying part of this 
equation, determining dynamics of the fields , can be separated out by averaging 

it according to Eq. (j3J). Under this operation, all terms linear in v±,pi give rise to o(l) 
contribution. For instance, 



(i) 



Of 



L 3 T 



d x v\ 



(i) 



to+T 



to 



^(5 



Lr 



T 



oil) 



(19) 



v 



1 If the gravitational field is not homogeneous, it is assumed to vary noticeably over distances larger than Lq. 
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in view of the estimates (|5]l. and the choice of L,T. The same argument applies to 
rfV)^, as well as to the second and fourth terms in the right hand side of Eq. (fTHj) . 
Furthermore, 



Vl (0) vWS 0) > = o 



according to the definition of V\ . Finally, contribution of the last term in the left hand side 
of (jl8j) also is o(l). Indeed, integrating by parts and taking into account Eq. (j!5|) . we have 



O^V) v?) = ^jdtjd*x („?>v) = ^fdtj(ds v?) v 



(i) 
1 5 



(20) 



where 5* is the surface of the cube V = {x : x G (xq,Xq + Ax), Axi = L}, ds being its 
element. Using Eqs. (JS)), the right hand side of Eq. (|2*Uj) is estimated as O(L /L). Hence, 



(«! 0) v) v{ 



o.)\ = Ll (h\ =oil] 



U \ L 



Thus, Eq. (fTHj) reduces upon averaging to the ordinary Euler equation for the functions 



(o) (o) 
*>o , Po 



A (0) , 

^- + (4°'V) 4 ' = -ivpf + G + 0(1) , (21) 

which proves the first part of the decoupling theorem. It is worth of mentioning that the 
large scale flow dynamics in the bulk turn out to be ideal at zeroth order in e. 

In connection with Eqs. (|15)) . (fTKj) the following circumstance should be emphasized. 
These equations describe bulk dynamics of the small scale parts of the flow variables at 
zeroth order in Lf/L Q , i.e., when the influence of the large scale flows on the flame cellular 
structure is completely neglected. This might seem to be in contradiction with the structure 
of Eq. (fTHjl . because it involves explicitly. However, the functions v^fa, t), pf (x,t), 
satisfying Eq. (fT7)|) in a given space-time region {x,t : x G (x ,x + Ax), Axi = L, 
t € (t ,t + T)}, can be written as 

v i ] ( x , t) = v[ 0) (x - v { 0) (x , t )t, t^j , pf ] (x,t)= pS 0) (x - v { 0) (x , t )t, , 

where , satisfy 

»~(°) 1 

^- + (v[°W) «<°> = -ivM 0) + ePrAv® . 
at V / q 

In other words, the role of in Eq. ()16|) is purely kinematical: it describes the large scale 
"drift" of the flame cellular structure. 
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C. Decoupling of jump conditions 



The proof of decoupling of the jump conditions is more complicated, since this is the place 
where the transport processes inside the flame front come into play. These conditions express 
the conservation of energy and momentum across the flame front. For freely propagating 
flames, and within the accuracy of the first order in the small front thickness, they were 
derived in the most general form in Ref . . To take into account the influence of gravity, it 
is sufficient to note that the bulk equations (jHj) can be rendered formally free by substituting 
p = p + g(Gx). However, gravity reappears through the jump conditions at the flame front. 
On the other hand, the influence of gravity on gas dynamics inside the flame front is small in 
comparison with the transport effects; their relative value is known to be given by the inverse 
Froude number Fr _1 = Lf ||</||/i7 f 2 = e. To the leading order in e, therefore, contribution of 
the gravitational field to the jump conditions is the same as in the case of a zero-thickness 
flame. 

For simplicity, we will consider two-dimensional (2D) case, assuming also that the Lewis 
number is equal to unity. No assumption is made concerning the incoming flow, except that 
its characteristic length L > L . Let the flame front position be described by an equation 
z = f(x,t), where the Cartesian coordinates (x,z) are scaled on Lq, and chosen so that 
z-axis is parallel to G. The x- and z-components of the flow velocity will be denoted by w 
and u, v = (w,u). We also introduce the unit vector r tangential to the flame front, and n 
orthogonal to it (pointing to the burnt matter). In components, 



Rewriting Eqs. (5.32)-(5.43) of jsj in this notation for the 2D case, and taking into account 
the contribution of the archimedean force to the pressure jump yields 




(22) 




(v+t) - («_r) 



) 



(23) 
(24) 



p+-p- 



P+-P- + (~qG z z) - {-gG z z)_ 



+ 




1 dfdN 
N dx dx 



) 



(25) 
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where 




G 



G 



Z ) 



and the subscripts "— " and "+" mean that the corresponding quantity is calculated for 
z = f(x, t) — and z = f(x, t) + 0, respectively. 

Finally, to complete the system of hydrodynamic equations and jump conditions, one 
needs an expression for the local burning rate. This expression, the so-called evolution 
equation, has the following form (Cf. Eq. (6.1) in Ref. P]) 



At this point, it is worth to make the following comment on the meaning of the above 

n 

asymptotic relations. Equations (|23j) - (|26|) were derived in Ref. |5J under assumption that 
the terms in the right hand sides of these equations, proportional to e, are small, which 
is only true if the gas flow is characterized by a length scale much larger than the flame 
front thickness. However, the rapidly developing LD-instability makes any smooth flame 
configuration highly corrugated within the time interval of the order As a result, the 
e-terms turn out to be of the order Lf/A c = 0(1). Similarly, account of the e 2 -corrections 
in the above equations would give rise to terms of the order Lf/X^ = 0(1) etc., questioning 
thereby validity of the small ^-expansion. However, it was mentioned in the Introduction 
that in practice, the cut-off wavelength A c is noticeably larger than the flame front thickness 
Lf. Thus, in the regime of fully developed LD-instability, the right hand sides of Eqs. (J23Tl - 
(|26j) are to be considered the leading order terms of the asymptotic expansion in powers 
of Lf/A c = 1/ci , rather than L{/Lq = e. On the contrary, e is the true parameter of the 
power expansions (|B|L (|7j). which determines the relative order of successive terms in these 
expansions. 

In order to extract from Eqs. jump conditions for the quantities Vq , , we 

need to introduce an auxiliary operation of averaging along the front. Let a quantity A be 
defined on the flame front, i.e., for {x, z, t : z = f(x,t)}. Given a point xq, choose Ax = 
Ax(xo) such that the front length £(t) between the points (xq, f(xo, t)) and (xo, f(xo+Ax, t)) 
satisfies 

A c < £(t) < L , £(f) = 0{L) 

for all t G (to, to+T). This is always possible since £(t) is of the order of distance between the 
two points. Then the average value of A over {x, z,t : x £ (x , x + Ax), t e (t , t + T),z = 




(26) 
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f(x, t)} is defined as 

t +T xo+Ax t +T to+T x +Ax 

(A) i = J dt J dl A, W= J dt £(t) = J dt J dl, (27) 

to XO to to Xo 

dl being the front line element, dl = Ndx. 

Using the operation introduced, the quantities r, n defined on the front, as well as the 
flame front position itself, can be decomposed into two parts corresponding to the scales Lq 
and A c , in a way analogous to Eqs. ©-(JZJ): 

t = t q + t 1 , (n x ,n z ) = (-t z ,t x ) , / = /o + /i, (t l ) i = 0, (A>i = 0, (28) 
ro,i = rg + 0(e) , /o = / (0) + 0(e) /i = 0(e). (29) 

To obtain jump conditions for the quantities expansions ©-(J7J) and (J2EJ should 

be inserted into Eqs. (J2"3|) - (j25|) . with the subsequent averaging of the latter along the flame 
front. 

Let us begin with the jump of the normal component of the gas velocity. In view of 
Eq. (jl3j) . one can introduce the stream function ip = ip(x, z, t) according to 

dtp dip 
M = ^' W = -^-z- (30) 

Using the operation of the bulk averaging (j3J), the function ip can be decomposed as 

V = V>o + ^i, (^i)=0. (31) 

It follows from Eqs. (jSJ) that 

1> = O(l), ^ = 0(e). (32) 
Substituting Eq. (J3T|) into Eqs. (jBD}, and averaging gives 



/d(iPo + ipi)\ diPo , 1 



t +T z +L 



Uo = \ 9x / = ^ + OT J dt J 

to zo 



n 



xo+L 



dx \ L 



and analogous equation for wq. Thus, up to o(l) terms, one has 



dipo dipo 
-o = ^, -o = -^- (33) 
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Inserting Eqs. (|50|l into Eq. (J23*j) . and using Eq. one finds 



t +T x +Ax 



<(w±w)>j 



2U 
1 

2H 



So 

t +T 



J-o 



dt ip { 



o± 



fa 



9-0 9-0 df 
dx dz dx ) ± 2U 



f +T x +Ax 

to XQ 



dx 



dip± 
dx 



(34) 



Since the variation of over space distances ~ L and time intervals ~ T is small, taking 
into account Eqs. (|29p. (jUS)), and neglecting 0(e) terms, one can write 



1 

2£f 



to+T 



dt ij}, 



'o± 



to 



xo+Ax 



xo 



2H 



fafatfao + Ax, /o(ar + Ax, t ),t ) - ipo±{xo, fo(x , t ),t )) 



T fdiPo + diP df 



22J \ dx dz dx 
Substituting this into Eq. (|23*|) gives 



Ax 



± 



T 
2U 



uo± ~ w o±~^: ) Ac • (35) 



u 0+ - u - - (w 0+ - 



dfo 
dx 



Note that 



and similarly, 



df /dx 
TV 



to+T x +Ax 



■ml dt 

U) XQ 



dx 



TAx 
2U 



dx dx 



l) + o(l). 



TAx df 
W dx 



(36) 



+ o(l), (37) 



TAx 

n zo = -gQ- + o(l) . 
Hence, Eq. (|37)|) can be rewritten as 

(v 0+ n ) - (v -n ) = (6 - 1) + o(l) 
or, with the same accuracy, 



v (0) n (0) 



l) + o(l) 



The inverse norm of no 



t +T 



~ y ^ £(t) 



l n o| 



to 



A.n/l+['^£ 



U 



<9x 



(38) 



(39) 



(40) 
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has a clear geometrical meaning. Namely, (H — 1) represents the relative increase of the 
flame front length due to its small scale wrinkling. 

Let us turn to the examination of the remaining jump conditions. Because of the e-terms 
in the right hand sides of Eqs. (J24j) . (J25|) . which involve highly nonlinear combinations of the 
flow variables and the function f(x,t), there seems to be a very little hope that the large 
scale parts of the flow variables eventually decouple from their small scale parts describing 
flame cellular structure. Nevertheless, they do, as will be shown presently. 

According to Eqs. (jH|)- (|TT|) and analogous estimates for the space-time derivatives of 
f(x,t), 

| = 0(1), f = 0(l), (41) 
^ oflY S-O^V i^-oflV (42) 



dt 2 \e J dxdt \e J dx 2 \e / 

the e-terms are 0(1). Let us show first that the average value of the right hand side of 
Eq. ()24j) along the flame front is actually o(l). Using the evolution equation (|26jl . expression 
in the parentheses in the right hand side of Eq. (J24|) can be rewritten as follows: 2 

b w _ + f'Du. + ±bf 
. /' , , (f'w-Y , CQV , ,, , n> 



u 



f> [u--N-f\ n 2 -1 
= v,_ + f'u_ + - + v ,_ W >_ + i J ~ + ^^ 

/ -\ N' 

+u>_(u--N-f) + - 

= v>- + ^f 1 - (u-f)' + A , (43) 

where the dot and the prime denote differentiation with respect to t and x, respectively. 
Hence, taking into account the estimates (JHJ), (ED), one has 

1 (a.. + + lay 



2 It was mentioned after Eq. that the e-terms in the jump conditions represent the leading order terms 
of the asymptotic expansion in powers of 1/ci . In transforming these terms, therefore, one can use the 
evolution equation with the e-term omitted. 
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dx [w- + f'u- 



t +T 

ir T + / dt 

to 



-u-f + 



u 2 + w 2 



xq+Ax 



xo 



o + off 



In view of the choice of T, £, the right hand side of Eq. 
0(Lf/£) = o(l). Thus, averaging of Eq. gives 



(44) 

turns out to be O(eL /£) = 



(M-M^o(i). 



To the leading order in e, this equation can be written as 



v 



(°) ..(0) 



OH 



V, 



0- 



.(0) 



.(o) 



0(1) 



-(0) 



T (o) , (0) 



(45) 



It was mentioned in the end of Sec. Ill ±31 that in the approximation considered, the large 
scale flows do not affect the flame cellular structure. In particular, the two directions along 
the flame front, r and — r, are left equivalent. This implies that up to 0(A c /£) terms, the 



value of ((uj+ — Vi_)t^)i must be invariant under the change 



.(0) 



. Hence, 



«2 - 



-(0) 



-(0) 



0. 



Thus, taking into account definition of Ti, we have from Eq. ()45|) 



(°) 



04 



0- 



.(0) 



oil). 



(46) 



Consider next the pressure jump, Eq. ()25p. The last term in the right hand side of this 
equation can be transformed as 3 

f'N' 



f + 2w-f + w 2 _f" + 2DN 



N 



(fl ATI \ f ATI 

]<i+w - N,+L w)- L w 



f + 2N + 2(u>_/)' + 2(w_N)' + (w 2 _f)' - 2w'_u_ + 



f'N' 
N 



f + 2N- (wtfj + 2w_u'_ + 



f'N' 
N 



(47) 



As before, the first three terms give rise to o(l) terms upon averaging. The remaining two 
terms, however, do not reduce to the full x- or t-derivatives. Their contribution is, therefore, 



See the footnote 2. 
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0(1), in general. Notice that the fourth term is quadratic in the gas velocity. This fact can 
be used to show that its contribution is independent of the functions t?o-,Po-- Indeed, one 
has, to the leading order in e, 

t +T xo+Ax 

J dt J dx («,£> + «,£>) («£> + uty' 



N , 2B 



to x 
t +T xo+Ax 



^ J dt J dx (wi?l + w?l)u?l' + 0(e) 

to XO 

to+T t +T xo+Ax 

e 



[ dt \wi?2u?l] X0+AX + ^- [ dt [ dx wflufl' + 0(e) 
J L J xo 2XT J J 



to to Xo 

t +T xo+Ax 
to Xo 

Thus, up to o(l) terms, the average value of ew_u'_/N turns out to be independent of the 
functions . The quantities wf}, w[- describe variations of the fuel velocity along the front 
cell in zero order approximation with respect to L{/Lq, i.e., when the influence of the large 
scale flow on the local flame structure is completely neglected. In particular, the value of 
(ew-u'_/N)i is independent of the coordinate x as well as of the time instant t . Denote 
this constant by a±/2 . 4 



4 At first sight, (ew-u'_/N)i depends on the choice of orientation of the coordinate axes, while the scalar 
pressure jump must be independent of this choice. It is easy to see, however, that within the accuracy 
of the above calculations, (ew-u'_/N)i is actually coordinate- invariant. In fact, under rotations of the 
coordinate system, w- , w_ transform as 

u>- — > cosy u>_ + sin ip w_ , 
m_ — > — sin ip W- + cos ip U- , 

where ip is the rotation angle. This transformation leaves (ew-u'_/N)i unchanged: 

(ew-u'_/N)i — ► ^ J dt J (cosy w_ + sin 93 u-)d(— sin 93 W- + cosip u_) 
= (ew-u'_/N),+o(^y 

Definition of a\ can be written also in an explicitly invariant form: 

/ w_u'_ — u_w'_ 
a ^ £ \ N 
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It remains to find the contribution of the last term in Eq. (|47|). We have 

t X 

The quantity (f'/N)' is nothing but the flame front curvature k, 

f" 

u — i 

k ~ ATS' 



Using the definitions (J22J), (|2Sjl . one can write 

= ( (^) ) = {T ' z)l = {T ' Z ° )l + {<l)l = T ' Z ° + ^ ' 
According to Eqs. (jUJ), (g2J), r^ = 0(1), = 0(l/e). Thus, to the leading order in e, 

e<*>, = e (rj?'^ = O(l) . 

Similarly to v± , describes geometry of a front cell neglecting the influence of the large 
scale flame structure on it. Hence, the value of (rj° }i is independent of the particular choice 
of the point xq on the flame front (and of the time instant to). Denoting this constant by 
hi, we thus obtain from Eq. (|25|) the following expression for the jump of at the flame 
front 

Pol ~ Po- = - d -^Gff ] - (9 - 1) - TTi + e[(9 - 1) - \n9]h + elnflai + o(l) , (49) 



where 7Ti is another constant defined by 



in, _ ,n,\ 



It is independent of xo, to on the same grounds as ai,ki. 

Finally, we have to average the evolution equation (|2fiJ) . As before, contribution of the £ 
term on the right hand side of this equation is o(l). Furthermore, 

xo+Ax 



t +T 

iV at / , 2B 



dx f 



X{) 



TAxdfo 
+ 0(1) = ^^ + 0(1 ^ 



to 



Using Eqs. (|3^l -<|55 jl . we find 

(«o_no) - n^/o = 1 + o(l) . (50) 
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Equations (J3HJ), (JIB)). P^j). and (joHj) constitute the proof of the second part of the de- 
coupling theorem. To make this more transparent, let us rewrite these equations in a more 
convenient notation. The jump conditions for the large scale parts of the flow variables read 5 

(tf+n) - (o_n) =(0-l>U + o(l), (51) 
(ff + t) - (o_t) = o(l) , (52) 

1 



P- = -^_^f + n + o(l), (53) 



where 



(0) / ,, -I \ , 

n ° f L ^ f=/ (o) , *=^(ff, ?=«r 



1 - || T (o)|| - ^k'k; ' 11 11 ~ 11 11 ~ ' 

p = pf, n = -(9-1) -7n + e[(0-l) -ln^fci+eln^!. 
The evolution equation takes the form 

p-fi)-i = U+o(l). (54) 

As we have seen, in zero order approximation with respect to e, the quantities fci, oti, 7Ti, 
and il entering these equations are independent of the coordinate xq and the time instant 
to- These constants can in principle be calculated provided that the exact small scale flame 
structure is known. However, as it follows from the above equations, this information is 
actually unnecessary, because dynamics of the large scale fields are independent of the 
specific values of these constants. Indeed, since the gas pressure enters dynamical equation 
only through its gradient, it is determined up to a constant. Therefore, the constants ki, «i, 
7Tx, are irrelevant. Unlike these, however, the constant il has a direct physical meaning. As is 
seen from Eq. ()54j) . il plays the role of the effective dimensionless velocity of the curved flame 
propagation. According to our choice of units, the gas velocity is scaled on the adiabatic 
velocity of a plane flame front C/f which is also used to define the units of pressure, time, 
and length [see the definition (jl2j) of L }. In analyzing the large scale flame dynamics, it is 
more natural to choose il • U{ , rather than U{ , as the velocity unit. Then il disappears from 
the jump conditions at the flame front, as well as from the flow equations (JEJ), (|2*T]) in the 



The vector quantities denoted by Gothic letters are designated with arrows. 
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bulk, which thus take the form of the equations governing propagation of a zero-thickness 
flame in an ideal fluid at constant speed iU7f with respect to the fuel, f being the flame 
front position, while n and t the normal and tangential unit vectors to f , respectively. The 
decoupling theorem is proved. 



III. NONLINEAR FLAME STABILIZATION IN GRAVITATIONAL FIELD 



The proved theorem considerably widens the scope of issues in flame dynamics accessible 
for analytical investigation. One of the most important consequences of the decoupling 
theorem is that unlike the local cellular dynamics, the large scale flame dynamics can be 
investigated in the framework of perturbation expansion with respect to the flame front 
slope, provided that the external field exerts a stabilizing influence on the flame. This fact 
will be illustrated below in the case of a flame propagating in an initially quiescent fluid in 
the direction of the gravitational field. As is well known (see, e.g., Ref. |ll|), in this case 
gravity plays the stabilizing role at the linear stage of development of the LD-instability. 
Our aim will be to determine the role of the nonlinear effects, and to explore the possibility 
of a full stabilization of the curved flame front by the gravitational field. 



We will follow the general met 
front position, developed in Ref. 



rod of deriving weakly nonlinear equations for the flame 



12 1 . It consists in bringing the system of hydrodynamic 



equations together with the jump conditions at the flame front to the so-called transverse 
representation in which dependence of all flow variables on the coordinate in the direction 
of flame propagation (z) is rendered purely parametric, and then reducing this system to 
a single equation for the front position. The calculation in the presence of gravity is very 
similar to that in the case of a freely propagating flame. Therefore, derivation of the equation 
will be only sketched below, referring the reader to the work |12| for more detail. 

As was mentioned in the end of the preceding section, it is natural to take il£7f as the 
velocity unit, redefining accordingly the units of length, time, and pressure to 

\\g\\ 

respectively. For simplicity, designation of the flow variables as well as space coordinates 
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and time will be left unchanged. Then the bulk equations read 6 



divfc = 0, (55) 

?+(Vv)o = -^ + <»> ® = ijr> 11^11 = !• ( 56 ) 

Since the fuel is assumed initially quiescent, the flow is potential upstream, and the general 
solution of Eqs. (|55p. (|56|) can be readily written down. In the reference frame of an initially 
plane flame front, 

+oo 

u = dk Ukex-p(\k\z + ikx) , (57) 

— oo 

to = H(u-l), (58) 
u + $(p + <5*) + |(u 2 + tt> 2 ) =0, (59) 
W = t) x , u = o 2 , (5 = -<S Z . 

Here H denotes the Hilbert operator defined as 

Hexp(ifcx) = i sign (A;) exp(ikx) , fe 7^ , (60) 
sign(A;) = |— . 

The LD-operator $ is related to H by $ = — H • d/dx. Equation (|59|) is nothing but the 
Bernoulli equation written in the transverse form. 

Because of the vorticity produced by the curved flame front, the flow of products of 
combustion is not potential. Nevertheless, the following transverse relation between the 
flow variables downstream can be obtained from Eqs. (}55|) . (}56|) at the second order of 
nonlinearity 

u-en'-$(ep + iu - 9 l + ^ +to(u' + ^+A =o, (ei) 

Equations (|5] | -()54| ) . (|58jl. ()59p. and (joT|) constitute the closed system of equations describing 
flame dynamics in the transverse representation. It can be reduced to the following equation 



For brevity, the o(l) symbols will be omitted in what follows. 
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for the function f 



+ l)f + 2#$f + 9(9 - l)f + (9 - l)0$f + (9- ^ — ^ ) $(f 2 



+^-^(5(f ) 2 + (e + 1) (ff + f Hf) + ^$ + (Hf) 2 ) 

+ (3£ - 1)$ (f Hf) + (TO - 1 + ~) f Hf - (Hf + <5f ) Hf = . (62) 

The linear terms in this equation reproduce the well-known equation 

(9 + l)f + 20$f + 9(9 - l)f + (9 - l)0$f = , 

from which it follows that at the linear stage of development of the LD-instability, the 
gravitational field plays the stabilizing role in the case of flame propagation in the direction 
of (25 ((25 = +1), and destabilizing in the opposite case ((25 = —1). To determine the role of 
the nonlinear terms, let us assume that there exists a stationary regime of flame propagation. 
It should be stressed that this assumption concerns only the the large scale front structure 
described by the function f. The local cellular structure does not need to be stationary Then 
Eq. (Jo^j) simplifies to 

- l)f + {9- l)<S<£f +(e- $(f') 2 + tS(f) 2 = . (63) 



2 J y " 9 

It is not difficult to see that the nonlinear term proportional to exerts a stabilizing influence 
on the flame if (25 = +1. Indeed, if we take 

f(x,t) ~ e at sm(kx), (64) 

Eq. (J62j) can be roughly considered as a "dispersion relation" for the increment o. As we 
see, the nonlinear term decreases a if & = +1, and vice versa. Therefore, Eq. (J63)) can 
only have solutions if <3 = +1, and the question of whether the flame can be stabilized 
by the gravitational field reduces to the question of existence of nontrivial solutions to this 
equation. 7 



7 The last term in the non-stationary equation 1)62(1 is also proportional to 0. Unlike the other two, it has 
a destabilizing effect on the flame front in the case (3 = +1. Indeed, using the definition of the Hilbert 
operator (|60|l . and substituting expression (|64|l into this term, one finds 

e ~ 1 ^U^2 / 



-f Hf = — cr|fc| cos^(fcc) < 0. 

v a 

This term, however, is irrelevant to the issue of existence of stationary configurations, since it contains f. 
Whether the gravitational field has an overall stabilizing effect, or not, depends on solvability of the 
stationary equation 
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Equation is a nonlinear integro-differential equation with respect to = f . To solve 
this equation, we first transform it as follows. Let us rewrite the nonlinear term $(f) 2 
iterating Eq. with respect to f", i.e., substituting 

f = -f$f+o(f 2 ). 



One has 

$m 2 = _H- , , 



$(f) 2 = -H^2! = _2H(f f ) = ^H(fdf) + O(f) . 



Using the well-known identity 

2R{^Rijj} = (Rip) 2 - i) 2 , 

we find 

Hf ) 2 = -j {(Hf ) 2 - (f) 2 } + 0(f) . 
Hence, within the accuracy of the second order, Eq. ()63|) takes the form 

f + «$f + /?(f) 2 -7($f) 2 = 0, (65) 

where 

® a ® fa i - 1 )^ * f a ( ~ V 



e ' ^ # 2 (#- 1) V 2 y ' ' 2 (e-i) V 2 

In connection with the transformation performed, it is worth to emphasize validity of the 



weak nonlinearity expansion when applied to t 



re investigation of the large scale flame 



dynamics. As was shown in detail in Refs. 0, 1C| . this expansion turns out to be self- 
contradictory in the case of a freely propagating stationary flame treated in the framework 
of the thin front model. This fact can be seen directly from Eq. (J63j) with (J5 = 0, in which 
case this equation reduces to the equality of two quantities of apparently different orders - f 
and $(f) 2 - 8 Only if — > 1 does the weak nonlinearity expansion of stationary flames make 
sense, since then f = 0(6 - 1), f = 0((6 - l) 2 ), so both terms in Eq. flH3|) are 0((6 - l) 3 ) 
quantities. On the contrary, in the presence of the gravitational field, f can be treated as 
the first order quantity when (9 — 1) is not small, because Of and G5$f are of the first order 



Perhaps, it is worth to stress once more that despite similarity of Eq. I|t)2l) with & = to that obtained in 
Ref. [12|, its meaning is completely different. In the notation of Sec. Iff Cl the latter equation determines 
the function /^ , while Eq. - the function . 



21 



in this case, and nothing prevents their sum from being formally a second order quantity. In 
fact, (0 — 1) must be finite in the latter case, since the second order term [9 — l)(3$f would 
be the only second order term in Eq. otherwise. As we will see below, solutions to this 
equation turn out to be unbounded for 9 — > 1. 

Turning back to Eq. (|65|). let us show first of all that its solutions, if any, are non-periodic. 
Notice that if f(x) is a periodic function, then [H(f — f)](x), where f is the mean value of f, 
is also periodic with the same period [substraction of f is necessary, because H is undefined 
on constants, see Eq. JSUJ)]. Let us integrate Eq. over period. The first two terms in 
this equation give rise to zero: 



J dx f" = f'| = 0, 

= - J dx (H(f-f)}'= H(f-f) 



On the other hand, using unitarity of the Hilbert operator, we find 

d*{/?(f) 2 -7(ftf) 2 } = / dx (/3 - >y)(f)* > . 

Of course, the absence of periodic solutions could be inferred already from Eq. (|5l5)). It is 
seen that such solutions are forbidden by the positive definite nonlinear term proportional 
to (5. 

Non-periodic solutions can be found in the form of the pole decomposition for the function 

<P = f 

2P l 

6 = a > , (66) 

^ x - x k 

k=l K 

where the value of the amplitude a, as well as position of P pares of the complex conjugate 
poles Xk, are to be determined by substituting this decomposition into Eq. (|B5|). Using the 
definition of Hilbert operator, one can show that 

TjA ■ \^ sl g n ( Im Xk) 

k=l 

It is not difficult to verify that Eq. is satisfied by (|6^j) provided that 

1 



P + l 
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and Xk satisfy the following system of algebraic equations 

ia{0 + l) + 2Y P Sign(Im %k) + 7 Sign(Im Xl) = , k = l,...,2P. (67) 



Evidently, this system is only consistent if = +1. Indeed, in the case of a pare of complex 
conjugate poles xi, x% = x{, one has, assuming Im x% > 0, 

g-7 (0-1) 2 
Xl "a(/5 + 7 )" 20 ' 

which is consistent with the assumed positivity of Im x\ if = +1. It is not difficult to 
show that the same is true in the general case of arbitrary number of poles. Consider the 
imaginary part of Eq. (|fi7j) corresponding to the pole uppermost in the complex plane, and 
take into account that (3 > 7. The fact that the system ()67|) is inconsistent for (25 = — 1 does 
not mean, of course, that the nonlinear stabilization is impossible in this case. Investigation 
of such a possibility requires account of higher order corrections. 
In the case of P = 1, one has, furthermore, 

f/= 1 ( 1 | 1 \ 2 (x-ReXi] 



(3 + 7 \x — xi x — x\J (3 + 7 (x — Re x\) 2 + (Im xi) 2 
and therefore, 

c 0(0-1). r, , 2 (0-i) 4 i 

f = In < [x - x ) H > , x = Re xi . 

The two-pole solutions for the cases 6 = 5, 10 and xq = are shown in Fig. 2. 

Equation (J65)) is derived in the scope of the power expansion with respect to the flame 
front slope <fi. In the case of the two-pole solution, |0| takes its maximal value 

K = T^l 

at the points x m = xq ± Im x%. We see that the developed weak nonlinearity expansion is 
valid if 9 is not too close to unity. For realistic values of the expansion coefficient [9 = 5 — 10) 
1. 

Next, consider the four-pole solution (P = 2). It has the form 

f = d SLz3_ i n { [( x _ Re Xl f + (i m Xl ) 2 ] [( x _ Re X2 )2 + ( Im X2 )2j 1 _ 
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FIG. 2: Two-pole solutions of Eq. (|65j) for 6 = 5 (full line), and 9 = 10 (dashed line). Flames 
propagate downwards. 

Assuming that Im x\,2 > 0, one has the following equations for the position of poles xi, X2, 

Ju 3 — \t 4 — ^° 2 

! ^ +7 » + 2 {^ + ( ^ 7) (^b; + ^)}^ ' 
ta( " +7 » + 2 {^ ± i + ( ^ 7) (^b + ^;)}^ - 

Separating the real and imaginary parts, and rearranging yields three equations for the four 
quantities Re X\^ , I m x i,2 

i P + 1 lm(x 2 -x 1 ) | 1 1_ = 0| (69) 

/3 — 7 l^i— X2I 2 Im X2 Im xi 

Re(xr - x 2 ) f . ^ + 7 |2 + , ^"^J = . (70) 
\ |a;i — X2\ pi — ^2!/ 

It follows from Eq. (J7U)) that Re X\ = Re £2 • This solution describes the "confluence" of 
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150 



100 



f(x) 




FIG. 3: Two-pole (full line) and four-pole (dashed line) solutions of Eq. for 8 = 8. 



poles. Then the remaining Eqs. (jUHj) . (JoTJj) give 



Im x lt2 = - ( 1 + 2 |— ^ 



1 ± 



2/3 



(# 2 - + 1) 1 ± 



20 



e 2 + i 



The two-pole and four-pole solutions are compared in Fig. 3 in the case of 9 = 8 and Xq = 0. 

The pole confluence is in fact a common property of the solutions (|66j). To see this, let us 
take the real part of the equation with k corresponding to the rightmost pole in the upper 
half-plane. We have 

2P 



(Re Xk — Re xi) 



i=i 
i^k 



W - xi\ 







In view of the choice of k, the left hand side is the sum of non-negative terms. It can be 
zero only if Re Xk — Re X\ = for all I. 

The question of which configuration is realized in the given conditions requires carrying 
out the stability analysis of various pole solutions, and can be solved, of course, only on 
the basis of the general non-stationary Eq. ([62)1 . According to the definition of f, such an 
analysis is to be performed with respect to perturbations with wavelengths A ~ L . 



25 



IV. DISCUSSION AND CONCLUSIONS 



The large scale flame dynamics are independent of its local cellular structure in zero 
order approximation with respect to the flame front thickness. This is the main result of 
the work, proved in Sec. |HI The local flame corrugation only affects the value of the normal 
velocity Uf changing it to iXUf, where il > 1 describes increase of the front length due to its 
wrinkling. In the scope of the thin front model, Uf plays the role of an external parameter 
specifying the characteristic velocity of the problem under consideration. Thus, the overall 
effect of the local flame structure on its large scale evolution amounts to a renormalization 
of this parameter. For flames of practical importance, the il-factor is about 1.3 — 1.5. In 
fact, it is itC/f, rather than Uf, which is more convenient to measure experimentally, since the 
measurement of Uf requires special facilities to suppress development of the LD-inst ability, 
such as those used in Ref. jl3|. 

The decoupling theorem allows one to avoid the difficult issues arising in investigating 
flame dynamics at length scales of the order A c , and to go directly to scales characterizing 
the problem in question. This is particularly important in numerical simulations of the 
flame dynamics. The computational grid should be chosen so as to well resolve the flame 
cellular structure, which leaves a little space for investigation of larger scales because of the 
limitations of computational facilities. 

The decoupling theorem also opens the way for analytical investigation of the large scale 
flame dynamics. As an example, the nonlinear development of the LD-instability in the 
presence of the gravitational field was considered in Sec. Illl| where a weakly nonlinear non- 
stationary equation for the flame front position was obtained [Eq. (J62)) ]. This equation 
admits stationary solutions in the case of flame propagation in the direction of the field, 
which means that the gravitational field has a stabilizing overall effect in this case. The 
resulting stationary flame configuration turns out to be essentially non-periodic, and repre- 
sents a symmetrical "hump" in the direction of the flame propagation, with slowly decreasing 
logarithmic "tails." A complete investigation of the non-stationary equation will be given 
elsewhere. 
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